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Abstract 

The pending question of the existence of three-neutron resonances near the physical 
energy region is reconsidered. Finite rank neutron-neutron forces are used in Faddeev 
equations, which are analytically continued into the unphysical energy sheet below the 
positive real energy axis. The trajectories of the three-neutron S-matrix poles in the states 
of total angular momenta and parity J n = | and J n = § are traced out as a function 
of artificial enhancement factors of the neutron-neutron forces. The final positions of the 
S'-matrix poles removing the artificial factors are found in all cases to be far away from 
the positive real energy axis, which provides a strong indication for the nonexistence of 
nearby three-neutron resonances. The pole trajectories close to the threshold E = are 
also predicted out of auxiliary generated three-neutron bound state energies using the 
Pade method and agree very well with the directly calculated ones. 
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I Introduction 



Controversial experimental hints on the existence of three-neutron (3n) resonances [1-4] and 
more recently even on four-neutron resonances [5] occurred again and again. No definite ex- 
perimental conclusion on 3n resonances seem to exist. On the theoretical side we are aware 
of various attempts to investigate states of 3n's but also here the situation is pretty unsettled. 
Obviously the optimal theoretical present day insight would be if the most modern neutron- 
neutron (nn) forces from the set of high precision NN potentials together with 3n forces of 
the 27r-exchange type (Tucson- Melbourne [6] or Urbana type [7]) would be used and the en- 
ergy eigenvalues of the 3n Schrodinger equation according to resonance boundary conditions 
would be determined in the unphysical energy sheet adjacent to the positive real energy axis. 
Unfortunately we can not provide such an insight right now but think we can improve at least 
on the existing studies we are aware of. 

One of the possibly first studies on that list is a solution of the 3n Faddeev equation based 
on a Yamaguchi rank 1 nn force in the state 1 So [8]. Only the 3n states of total angular 
momenta J w = 1/2" and 3/2" (degenerate in this case) have been studied. The nn force was 
artificially enhanced such that two and three neutrons were bound. Then the enhancement was 
reduced which has the consequence that both bound state energies move towards E — 0. For 
the forces used it happened that the 3n binding energy moves faster than the 2n binding energy 
and thus hits the dineutron-neutron threshold to the left of E — 0. The trajectory of the 3n 
bound state energy continues then into the unphysical sheet below the dineutron-neutron cut. 
In this manner a 3n resonance occurs, below the 3n break-up threshold at E — 0. Further 
decreasing the enhancement factor the 3n resonance trajectory moves up again exactly towards 
E = where it meets with the In bound state energy. Both energies coincide exactly at E = 
corresponding still to an enhancement factor larger than 1. In further decreasing the potential 
strength till one reaches the physical value 1 the vanishing 2n bound state energy turns into the 
well established 2n virtual state energy in the second sheet on the negative real energy axis and 
the 3n resonance disappears in a unphysical sheet which is different from the sheet adjacent 
to the positive real axis. In other words the complex resonance energy of that 3n resonance 
for J w = 1/2" and 3/2" has no positive real part and therefore can not be felt for positive 
three-neutron energies. We have to emphasise that this refers to a nn force acting only in the 
state 1 S'o and it will change if additional nn force components will be added as will be shown 
in this paper. 

A very much related study has been performed in [9] with essentially the same result as in [8]. 
Using just bound state techniques [10] the necessary enhancement factors on the nn forces in 
S- and P-waves have been determined to bind 3 neutrons near zero energy. Based on the Reid 
potential [11] enhancement factors of the order of 4 have been found which make low lying 3n 
resonances quite unlikely. 

Another theoretical investigation [12] we are aware of for the state J w = 1/2" has been per- 
formed on the basis of a hyperspherical harmonic expansion. Local forces like the ones of 
Mafliet Tjon have been used and this only in the nn state 1 S . The expansion was truncated 
and the zeros of the Jost function were determined. It resulted a resonance energy around 
E = —4.9 — i 6.9 MeV in an unphysical sheet. The authors found a strong sensitivity to the 
choice of the potential. 

A further investigation [13] relied on a variational ansatz and used the complex scaling method 
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which allows to treat a resonance problem like a bound state one. The Minnesota effective 
nn force [14] together with tensor forces [15] have been applied. The authors find as the only 
candidate for a 3n resonance the state J 71 " = 3/2 + and a prediction for its energy of E = 14 — % 13 
MeV. 

Finally in a more recent work [16] the Faddeev equation in configuration space has been solved 
for the states J n = l/2~ and 3/2 ± using the more realistic forces by Gogny et al. [17] and 
the Reid93 [18] potential in the 2n states 1 5'o, 3 Po and 3 P 2 - 3 F 2 . Choosing proper boundary 
conditions according to 3n resonances the complex energy eigenvalues are determined starting 
again from artificially enhanced forces and reducing their strengths gradually. Unfortunately 
with increasing negative imaginary parts of the complex resonance energies numerical instabil- 
ities occurred and the trajectories could not be followed up until the physical values for the 
enhancement factors have been reached. 

Now we would like to improve on all that by solving the Faddeev equations with low rank nn 
forces in all relevant partial wave states exactly and determining the final resonance positions 
for the actual force strength as prescribed by nn phase shifts (assumed to be identical to the 
strong pp phase shifts). In the light of our results we shall also comment on the previous 
findings. 

This paper is organised as follows. In section II we briefly review the set of coupled Faddeev 
eigenvalue equations for a 3n system based on finite rank forces and mention the steps needed 
for an analytical continuation into the unphysical energy sheet adjacent to the positive real 
axis. The dynamical nn force input and the resulting 3n resonance trajectories for the states 
J 71 " = \j2 ± and 3/2 ± are presented in section III. In section IV we show how the part of the 
3n resonance trajectory close to the threshold E = can be predicted with the help of Pade 
approximant's from a set of 3n bound state energies. This serves as a test for our results 
obtained through the analytical continuation of the Faddeev equation. Though the final 3n 
resonance positions could not be reached with that method, we think this illustration should 
be of interest since it works beautifully for resonances close to the first threshold and requires 
only bound state techniques. We summarise in section V. 



II Formalism 

The Faddeev equation for three identical nucleons has the well known form [19] 

ij = G tP^j, (1) 

where Gq is the free 3N propagator, t the NN t-matrix and P the sum of a cyclical and anti- 
cyclical permutation of 3 objects. We regard that equation in a momentum space representation 
and partial wave decomposed. It results a set of coupled equations in two variables, p and q, 
which are the magnitudes of two relative momenta. In the present investigation we shall use t- 
matrices of finite rank. The set of two-dimensional equations changes then into a set of coupled 
one- dimensional ones. They have the well known form 

Q 3 a{<l) = I d Q'q f2 J2Yl I dx^-^-G aa ,(q,qf,x) 

J ° ocl km J - 1 n l 

£ ^ £ _ ^ j£ 2) tL(z>)<C(<1<) (2) 
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and 6 is connected to ip via 

i> a (p,q) = — ^E^riii^f?). (3) 

^ M AM 3,k 

The underlying two-body t-matrices are 

t{p, p', z) = Y, 9i(p) Ti,j(z) 9j(pf) (4) 

The corresponding form factors g{p) depend on the obvious two-body quantum numbers a = 
(Isj) . The two-body subsystem energy z = E — |^ is given in terms of the total energy E and 
the kinetic energy of relative motion of the third particle. Further a stands for the string of 
three-body quantum numbers 

« = ^)j{\\)I{jI)J{t l -)T (5) 

For all those by now standard notations we refer to [19]. Finally G aa ,(q,q/,x) are geometrical 
coefficients arising from the partial wave representation of the permutation operator P and 7Ti 
and 7r 2 are shifted arguments related to P and depend on q, ql and x, see [19] for details. 
We are interested in 3n bound state and resonance energies, the latter ones emerging from 
bound state energies when the artificially increased interaction strength is reduced. Since two 
neutrons are not bound in nature the only unphysical sheet of interest is the one which can 
be reached through the cut along < E < oo. It will be denoted in the following simply 
as the unphysical sheet. While for E < there is no singular denominator, for E > the 
free propagator in Eq. (2) is not defined. The well established manner to analytically continue 
the coupled set of equations in the energy into the unphysical sheet is to deform the path of 
integration into the lower half plane. This has been described for instance in [8] and [20]. The 
singularities which are avoided in this manner arise from the free propagator and the form 
factors. The r-function is not singular for energies z = E — j n the unphysical sheet under 
investigation. A typical path in the complex g-plane (for both q and ql) together with domains 
of singularity arising from the free propagator for the example E = 10.0 — % 10.0 MeV is shown 
in Fig. 1. 

The strategy now as pointed out in the introduction is to first enhance the nn forces such that 
a 3n bound state exists. At the same time also two-neutron bound states might exist. Then by 
suitably chosen enhancement factors one can achieve that the two-body bound state energies 
coincide and move together towards E — 0, when the enhancements are reduced. In all cases 
we studied the 3n bound state energies stayed always to the left of the two-body bound state 
energies until they reached E — 0. This was not the case in [8], where for the degenerate 
states l/2~ and 3/2" only the 1 S'o force had been kept. Here we include on top P- and D- 
wave forces and such a coincidence of 3n and 2n bound state energies does not occur. Once 
the two-neutron bound state energies have disappeared we still have a 3n bound state. Then 
reducing the enhancement factors further the 3n bound state energy will approach E — and 
then the energy eigenvalue E will dive into the lower half plane of the unphysical sheet. There 
we shall follow its trajectory until the physical strengths of the nn forces are reached. This 
defines the final positions of the 3n resonances. It is well known and easily seen that what we 
call 3n resonance energies are poles of the 3n ^-matrix. 
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Figure 1: The regions of singularities from the free propagator for the complex energy E = 
10.0 — i 10 MeV and a path of integration C. 

Ill Results 

We use nn forces of rank 2. They are given in Appendix A. They describe pp phase shifts 
(without Coulomb force) reasonably well as documented in the Appendix. We restrict our 
study to the states 1 S , 3 P , 3 P2 and 1 D 2 . The coupling between 3 P 2 and 3 F 2 is neglected. 
The nn force in the state 3 Pi is repulsive and therefore we did not include it in the main 
investigation, since it has to be expected, that it will move the 3n resonance position further 
away from the positive real energy axis. At the end we shall consider its effect and show that 
its influence has indeed repulsive character but only a marginal one and could therefore indeed 
be neglected. 

We shall investigate now the 4 possible states J 71 " = l/2 ± and J 71 " = 3/2 ± in turn. 
3.1 The state 3/2 

We use up to 10 three-neutron channels which are displayed in Table 1. To reach a 3n bound 
state it is sufficient to enhance just the 3 P 2 nn force and to keep all the other force components 
at their physical values. As an example for an enhancement factor we quote \-ap 2 = 3.5 where 
the 3 neutrons in a 10 channel calculation are bound at —7.79 MeV. The threshold energy E = 
is reached for \sp 2 = 3.23. For the sake of future comparisons we provide a few intermediate 
3n resonance positions in Table 2. The 3n resonance trajectory is shown in Fig. 2 for a 10 
channel calculation. The trajectories for a smaller number of channels are very similar. The 
final resonance positions for an increasing number of channels in the order given in Table 1 
are displayed in Table 3. Adding more nn force components will certainly not change the final 
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Table 1: The partial wave quantum number for the three neutron state J n = | 



nn state 


/ 


s 


j 


A 


/ 


'So 











1 


3/2 


3 P2 


1 


1 


2 





1/2 


3 P 


1 


1 





2 


3/2 


3 P 2 


1 


1 


2 


2 


3/2 


3 P 2 


1 


1 


2 


2 


5/2 


^2 


2 





2 


1 


1/2 


l D 2 


2 





2 


1 


3/2 


l D 2 


2 





2 


3 


5/2 




1 


1 


2 


4 


7/2 




2 





2 


3 


7/2 















position in a significant manner in the sense that the position would come somewhere close to 
the positive real energy axis. If one adds the 3 Pi nn force, which is of repulsive nature, we have 
13 channels and the final resonance position shifts by —1.79 — i 2.87 MeV, which is marginal. 
In Ref.[16] which can be considered to be the most realistic approach towards three-neutron 
resonances carried out so far the rather realistic nn forces by Gogny [17] and Reid93 [18] have 
been used. In contrast to our simplified forces they also include the tensor force induced coupling 
3 P 2 - 3 F 2 . For the case of the Gogny potential the enhancement factor for the 3 P 2 — 3 F 2 nn 
force component was about 4.4 to have 3n's bound at zero energy, while for our simplified force 
this was the smaller value of 3.5. In [16] reducing A3p 2 _3^ 2 to 3.7 yields a resonance energy of 
4.95 — % 3.75 MeV, while at this value our 3n system is still bound. Thus one might conjecture 
that in the more realistic case the 3n resonance would be even further shifted away from the 
positive real energy axis. On the other hand the motion of our resonance energy for a change 
of AA = 4.4 — 3.7 = 0.7 from the threshold E = (corresponding to an enhancement factor 
2.8) is 5.58 — i 4.04 MeV, which is slightly larger. Therefore a direct quantitative comparison 
is not possible, though qualitatively the results are similar. This is also true for the case of 
Reid93, also considered in [16], where the enhancement factor \z P2 _z F2 = 3.25 yields a resonance 
position of 5.30 — i 3.53 MeV. 



3.2 The state 3/2+ 

Again we use up to 10 channels, which are displayed in Table 4. In order to have a 3n bound 
state we had to enhance the 3 Pq and 3 P 2 nn forces while the 1 S'o and 1 D 2 forces could be kept 
at their physical values. For the enhancement factors \sp 2 = 3.5 and \3 Po = 5.0 2n bound 
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Table 2: Intermediate 3n resonance positions together with the enhancement factor A3p 2 for a 10 
channel calculation and J w = | . 



A 3 p 2 


EIT [MeV] 


3.22 


0.20 -i 0.003 


3.20 


0.57 -i 0.03 


3.00 


3.51 - i 1.45 


2.00 


5.40 -i 19.17 


1.50 


-1.95 -i 28.87 


1.00 


-12.13-i 37.96 



T 



i , i , i , i , i 

-10 -5 5 10 

Re E [MeV] 

Figure 2: The resonance pole trajectory for the state 3/2~. 
Table 3: The final 3n resonance positions for the state J 71 " = | including different number of channels. 



Number of channels 


Final resonance positions 


2 


-12.81-i 33.22 


3 


-13.14-i 33.21 


4 


-13.11-i 33.17 


5 


-13.33-i 32.91 


6 


-14.28-i 33.37 


7 


-12.43-i 37.65 


8 


-12.31-z 37.74 


9 


-12.31-z 37.74 


10 


-12.13-2 37.96 



-20 



-35 
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Table 4: The partial wave quantum number for the three neutron state J n = | 



nn state 


/ 


s 


7 


A 


/ 


'So 











2 


3/2 


3 Po 


1 


1 





1 


3/2 


3 P 2 


1 


1 


2 


1 


1/2 


3 P 2 


1 


1 


2 


1 


3/2 


l D 2 


2 





2 





1/2 




1 


1 


2 


3 


5/2 


3 P 2 


1 


1 


2 


3 


7/2 


^2 


2 





2 


2 


3/2 


^2 


2 





2 


2 


5/2 


^2 


2 





2 


4 


7/2 















states for 3 P 2 and 3 Po just disappear and the 3 neutrons are bound at —6.55 MeV. Then we 
reduce both enhancement factors. The way this is done is in principle arbitrary. Again for 
future work we provide some intermediate 3n resonance positions along the trajectory in the 
unphysical sheet in Table 5. The corresponding trajectory is shown in Fig. 3. The dependence 
of the final position on the number of channels is displayed in Table 6. Like for 3/2 ~~ the final 
^-matrix pole position is far away from the positive real energy axis. Again the effect of 3 Pi 
nn force is marginal (a shift of —0.80 — i 0.42 MeV). 

Now we compare to the work in Ref. [16]. The enhancement factors for the Gogny potential for 
the 3 P ari d 3 P 2 — 3 F 2 nn force components are about 7.7 and 4.9 to generate nn bound states 
at zero energy. This is larger in both cases than for our simplified potential. Correspondingly 
according to Ref. [16] for \3 P{) = 5.0 and A3p 2 _3 Fa = 4.30 one ends up with a resonance at 
5.74 — i 1.53 MeV, while for such enhancement factors our nn forces still bind the 3n system. 
Thus also here one is tempted to conjecture that the more realistic Gogny potential would 
lead to a final resonance position, which is farther away from the real axis than what we find. 
For the Reid93 potential the corresponding values are 5.93 — i 1.55 MeV for A3 Po = 4.6 and 
\ip 2 -AF2 = 3.50. 

3.3 The state 1/2 

The partial wave quantum numbers for that state are shown in Table 7. In this case there are 
only 5 channels. It was necessary to enhance the 3 P 2 and 1 D 2 nn forces to reach a 3n bound 
state. For A3p 2 = 3.5 and Ai£> 2 = 6.90 the In bound states just disappear and the 3n system is 
bound with —1.90 MeV. Again we display some intermediate resonance positions in Table 8. 
The final S'-matrix pole position is significantly further away from the positive real energy axis 
than for the states 3/2 ± . The resonance trajectory is shown in Fig. 4. If one include 3 Pi, one 
has 7 channels and the resonance position shifts by —4.65 — i 0.53 MeV which is insignificant 
in relation to the value without 3 Pi. 

This case poses even more difficulties to compare to Ref. [16], since there 3 Pq and 3 P 2 — 3 F 2 have 
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Table 5: Intermediate 3n resonance positions together with the enhancement factor A3 P and X^ p for 
a 10 channel calculation and J w = § . 



Asp 2 


A 3Po 


Eg- [MeV] 


3.20 


4.70 


0.60-i 0.004 


3.00 


4.50 


3.21-i 0.96 


2.00 


3.50 


6.41-i 14.75 


1.80 


3.30 


5.12-i 18.18 


1.00 


2.50 


-3.53-i 29.96 


1.00 


1.00 


-10.30-i 33.20 



Table 6: The final 3n resonance positions for the state J n = | including different number of channels. 



Number of channels 


Final resonance positions 


3 


-9.82-i 58.37 


4 


-10.06-i 58.31 


5 


-10.24-i 33.16 


6 


-10.25-i 33.17 


7 


-10.61-i 33.37 


8 


-10.46-i 33.29 


9 


-10.35-i 33.23 


10 


-10.30-i 33.20 



been enhanced while here we found it more natural to enhance 3 P 2 and 1 D 2 (the latter force 
component was not included in Ref. [16]). For the enhancements A3p = 5.5 and A3p 2 _3p 2 = 2.75 
in [16] the 3n resonance turned out to be 2.82 — % 2.82 MeV. Apparently these factors are still 
quite large that also here one cannot expect a nearby resonance. 



3.4 The state 1/2+ 

Here we used 6 channels displayed in Table 9. Again we need to enhance both the 3 P 2 and 
1 Z?2 nn forces by 3.5 and 6.90, respectively, to get a 3n bound state at —7.38 MeV. Some 
intermediate resonance positions are shown in Table 10. The corresponding trajectory is shown 
in Fig. 5. The inclusion of 3 Pi leads to 8 channels and makes a tiny shift to the final resonance 
position of —0.36 — % 0.02 MeV. In Ref. [16] this state has not been considered. 
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-10 -5 5 

Re E [MeV] 

Figure 3: The resonance pole trajectory for the state 3/2 + . 
Table 7: The partial wave quantum number for the three neutron state J n 



nn state 


I 


s 


3 


A 


I 


% 











1 


1/2 


3 Po 


1 


1 








1/2 




1 


1 


2 


2 


3/2 


3 P2 


1 


1 


2 


2 


5/2 


l D 2 


2 





2 


1 


3/2 



Table 8: Intermediate 3n resonance positions together with the enhancement factor and 
for a 5 channel calculation and J n = \ . 



Xsp 2 


Ai Da 


EIT [MeV] 


3.38 


6.88 


0.09-i 0.0002 


3.00 


6.50 


7.32-z 3.49 


2.00 


5.50 


16.44-i 23.16 


1.50 


5.00 


14.51-i 32.15 


1.00 


4.00 


6.10 -i 42.27 


1.00 


3.00 


-2.25 -i 46.45 


1.00 


2.00 


-11.36-i 47.65 


1.00 


1.00 


-20.38-i 45.27 



10 



J 

-20 -10 10 20 

Re E [MeV] 

Figure 4: The resonance pole trajectory for the state 1/2" . 
Table 9: The partial wave quantum number for the three neutron state J n = \ 



nn state 


I 


s 


3 


A 


J 


% 














1/2 


3 Po 


1 


1 





1 


1/2 


3 P 2 


1 


1 


2 


1 


3/2 


3 P 2 


1 


1 


2 


3 


5/2 


l D 2 


2 





2 


2 


3/2 


l D 2 


2 





2 


2 


5/2 



Table 10: Intermediate 3n resonance positions together with the enhancement factor A3p 2 and Ai£> 2 
for a 6 channel calculation and J 71 = i + . 



Aap 2 


^D 2 


EIT [MeV] 


3.43 


6.43 


0.64 -i 0.002 


3.00 


6.00 


12.90-i 4.46 


2.00 


5.00 


16.91-i 22.63 


1.00 


4.00 


11.94-2 34.46 


1.00 


3.00 


3.97 -i 41.41 


1.00 


2.00 


-4.55 -i 45.51 


1.00 


1.00 


-14.29-i 48.34 
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Figure 5: The resonance pole trajectory for the state l/2 + . 

IV Resonance Prediction via Pade 

The interesting idea to predict low lying resonances with the help of Pade approximants from 
a sequence of (auxiliary) bound state energies has been proposed in [22] and exemplified for 
instance in [23]. We also applied it to the system of 3n's. Assume, the 3n binding energy is 
given in terms of a power series in an enhancement factor A as 

E(X) = E(X ) + (A - A„)*ai + (A - A ) €+1 a 2 + +(A - A ) ?+w - 1 a m (6) 

where £ has still to be determined. In searching for resonances close to the threshold E = it 
is natural to choose A such that E(\ ) = 0. For two particles in this case it is known [22] 
that £ = 1 for a S- wave and £ = \ for P-waves and higher ones. For a three-body system 
and the three-body break-up threshold we are not aware of an analytical insight into the value 
£. We determined £ numerically by approaching E = from both sides, A greater than Ao 
and A smaller than Ao- We found £ = | within our numerical accuracy for all studied cases 
J n = l/2 ± and J 71 " = 3/2 ± . Choosing m 3n bound state energies close to E = one can 
determine a\, Ch2, ■ ■ ■ a m . To do that we solved the resulting m equations by the method 
of singular value decomposition [24]. The power series (6) can then be rewritten in form of 
Pade approximants. We used a continued fraction expansion, which is equivalent to certain 
orders of approximants (see for instance [19]). Choosing for instance m = 6 one can reach the 
Pade approximant of the order [3,3]. This turned out to be sufficient to predict the low lying 
resonances displayed in Tables 12,14,16 and 18. In addition we show the bound state energies 
in Tables 11,13,15 and 17 out of which the a m have been determined. There we also include the 
resonance positions determined directly from the analytically continued Faddeev equation. In 
the cases J n = | + and ~ we changed two enhancement factors at the same time, but exactly 
by the same amount AA. In this case AA is the expansion parameter in Eq. (^). The a m 's as 
well as the resonance energies are then determined either through the enhancement factors in 
one or the other state. The resulting resonance positions are exactly the same. The last rows 
in Tables 13, 15 and 17 show the two Ao-values to many digits required to reach a significantly 
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Table 11: The J 71 " = | bound state energies used as input for Pade [3,3]. The last row shows Ao and 
the resulting energy. 



\$p 2 


fpSn nVTpVl 
B l-*-*-'-^ v i 


3.41 


-4.79366 


3.36 


-3.28446 


3.31 


-1.89909 


3.26 


-0.65666 


3.25 


-0.42848 


3.24 


-0.20848 


3.2300873 


-l.lxlO" 8 



Table 12: The J n = § 3n resonance positions determined directly from the Faddeev equation in 
comparison to the Pade [3,3] prediction. 



Asp 2 


E^ s [Me] via Faddeev 


Pade[3,3] [MeV] 


3.22 


0.20-2 0.003 


0.20-i 0.003 


3.21 


0.39-i 0.014 


0.39-i 0.012 


3.20 


0.57-2 0.03 


0.57-i 0.03 


3.19 


0.74-2 0.06 


0.74-i 0.05 


3.18 


0.91-2 0.09 


0.90-i 0.08 


3.17 


1.08-2 0.13 


1.06-i 0.12 


3.16 


1.24-2 0.18 


1.21-i 0.16 


3.15 


1.40-i 0.23 


1.35-i 0.22 


3.14 


1.56-2 0.28 


1.50-i 0.28 


3.13 


1.71-i 0.34 


1.63-2 0.35 



small 3n binding energy. In all cases we see a very good agreement of the low lying resonances 
using the Pade method with the directly calculated ones. It is important to note that the 3n 
bound state energies which serve as input have to be known to 5 digits in order to achieve the 
quoted Pade results. 
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Table 13: The J n = | bound state energies used as input for Pade [3,3]. For the last row see text. 



Xap 2 


A 3Po 


E\ n [MeV] 


3.41 


4.91 


-4.04286 


3.36 


4.86 


-2.76896 


3.31 


4.81 


-1.58830 


3.26 


4.76 


-0.51213 


3.25 


4.75 


-0.31107 


3.24 


4.74 


-0.11532 


3.2339687436 


4.7339687436 


-1.4xl0" 9 



Table 14: The J n = § 3n resonance positions determined directly from the Faddeev equation in 
comparison to the Pade [3,3] prediction. 



Asp 2 


■Wo 


E% s [MeV] 


Pade [3,3] [MeV] 


3.22 


4.72 


0.26-2 0.0001 


0.26-i 0.0006 


3.21 


4.71 


0.43-i 0.001 


0.43-i 0.003 


3.20 


4.70 


0.60-i 0.004 


0.61-i 0.007 


3.19 


4.69 


0.76-2 0.012 


0.77-i 0.013 


3.18 


4.68 


0.92-2 0.024 


0.93-i 0.022 


3.17 


4.67 


1.06-2 0.04 


1.08-i 0.03 


3.16 


4.66 


1.21-2 0.06 


1.23-i 0.05 


3.15 


4.65 


1.35-2 0.09 


1.37-i 0.07 


3.14 


4.64 


1.49-2 0.13 


1.50-i 0.09 


3.13 


4.63 


1.62-i 0.17 


1.63-2 0.11 


3.12 


4.62 


1.76-2 0.21 


1.76-i 0.14 



Table 15: The J n = \ bound state energies used as input for Pade [3,3]. For the last row see text. 



Aap 2 


Al£> 2 


E 3 B n [MeV] 


3.49 


6.89 


-1.58303 


3.48 


6.88 


-1.27353 


3.47 


6.87 


-0.96991 


3.46 


6.86 


-0.67262 


3.45 


6.85 


-0.38224 


3.44 


6.84 


-0.09983 


3.436380542 


6.836380542 


-1.8xl0 -7 
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Table 16: The J n = \ 3n resonance positions determined directly from the Faddeev equation in 
comparison to the Pade [3,3] prediction. 



A 3 p 2 


Al£> 2 


E'% s [MeV] via Faddeev 


Pade [3,3] [MeV] 


3.43 


6.83 


0.17-i 0.001 


0.17-i 0.001 


3.42 


6.82 


0.43-i 0.009 


0.43-i 0.01 


3.41 


6.81 


0.69-i 0.02 


0.69-i 0.03 


3.40 


6.80 


0.94-i 0.05 


0.94-i 0.06 


3.39 


6.79 


1.19-z 0.07 


1.19-i 0.09 


3.38 


6.78 


1.43-i 0.10 


1.44-i 0.13 


3.37 


6.77 


1.67-i 0.14 


1.69-i 0.18 


3.36 


6.76 


1.91-i 0.18 


1.94-i 0.23 



Table 17: The J n = \ bound state energies used as input for Pade [3,3]. For the last row see text. 



A3p 2 


^D 2 


E 3 B n [MeV] 


3.50 


6.90 


-7.37894 


3.45 


6.85 


-5.02695 


3.41 


6.81 


-3.22567 


3.36 


6.76 


-1.08684 


3.35 


6.75 


-0.67596 


3.34 


6.74 


-0.27141 


3.333194193 


6.733194193 


-2.0xl0~ 7 
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Table 18: The J n = ^ 3n resonance positions determined directly from Faddeev 
equation in comparison to Pade [3,3] prediction. 



Asp 2 


Al£) 2 


E^s [MeV] via Faddeev 


Pade [3,3] [MeV] 


3.33 


6.73 


0.13-i 0.0000 


0.13-i 0.0002 


3.32 


6.72 


0.52-i 0.0007 


0.52-i 0.003 


3.31 


6.71 


0.90-i 0.004 


0.90-i 0.007 


3.29 


6.70 


1.27-i 0.009 


1.28-i 0.011 


3.28 


6.69 


1.64-i 0.02 


1.66-i 0.01 


3.27 


6.68 


2.00-i 0.03 


2.02-i 0.02 


3.26 


6.67 


2.35-i 0.05 


2.38-i 0.02 


3.25 


6.66 


2.70-i 0.07 


2.73-i 0.03 


3.35 


6.65 


3.04-i 0.09 


3.07-i 0.03 



V Summary 

We analytically continued the Faddeev equation for three neutrons into the unphysical sheet of 
the complex energy plane adjacent to the positive real energy axis. The nn forces are of finite 
rank and were kept in the partial waves 1 S'o, 3 Pq, 3 Pi, 3 P2, 1 P > 2- This leads up to 13 channel 
calculations. The 3 Pi nn force is repulsive and its effect turned out to be rather small in relation 
to the combined effects of all the other forces. The analytical continuation is performed via a 
contour deformation of the path of integration. The strategy is to artificially enhance the nn 
forces such that a 3n bound state is generated. Then by reducing that enhancement one follows 
the energy eigenvalue of the homogeneous 3n Faddeev equation which traces out a trajectory. 
After having reached the 3n break-up threshold E — that trajectory dives into the unphysical 
energy sheet. The final positions for the physical strength turned out for the four cases studied, 
J n = l/2 ± and J n = 3/2 ± , to be far away from the positive real energy axis. Though the nn 
forces employed are not of the quality of the modern high precision NN potentials it appears 
unlikely that those modern forces would lead to totally different results. Thus we conjecture 
that 3n resonances close to the physical region will not exist. This conjecture is also supported 
by the findings in [16]. 

These results gained through the analytically continued Faddeev equation have been tested 
numerically using Pade approximants. In this method one generates a set of auxiliary 3n 
bound state energies close to E — for different enhancement factors A. This determines the 
first few coefficients in a power series expansion of E(X). This again leads to the coefficients of 
Pade approximants which can be used to predict low lying 3n resonances. The resulting values 
agreed very well with the ones determined via the Faddeev equation. 
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Table 19: The parameters of our nn potential. The units of these parameters are fm 1 for 
MeV fm" 2 ^ 1 for A« and the % are dimensionless. 







3p 

-TO 


:i Pi 


:i P2 


L D 2 


011 


0.8131678 


0.8322894 


0.9713441 


1.977127 


2.522169 


Pl2 


1.288463 


1.262785 


2.180297 


3.147108 


1.651999 


021 


7.496476 


2.645783 








P22 


1.661389 










7i 


2.698168 


6.397468 


31.97981 


5.242138 


-0.2668325 


72 


0.3270664 










Ai 


-17.87098 


-12.10082 


59.89249 


-1477.246 


-918408.7 


A 2 


82710.93 


277071.8 









Table 20: The nn effective range parameter r and the scattering length a in units of fm for the 
1 So-state in comparison with experimental results. 





r 


a 


our potential 
exp. data [21] 


2.85 
2.83±0.11 


-19.08 
-18.6±0.5 
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Appendix A 

The finite rank potential we use has the following form [25] 

V(p,pf) = j29i(p)\9M (1) 
i=i 

where 

p i p l+2 

3M = ~2 7 WWi + (r) 2 7 W ( 2 ) 



(p 2 + P 2 n) l+1 (p 2 + f3hY 
p l+2 p l+i 

(p 2 + Piy +2 + (p 2 + Pi 2 y 



9^(P) („2 i fl2 V+2 + ( n 2 i R2\l+3^ 2 ^) 



The /-dependent parameters of the potential form factors are displayed in Table 19. The 
resulting phase shifts of the nn system for the partial waves 1 S , 3 P , 3 Pi, 3 P 2 and 1 D 2 are 
compared to pp phase shift values (without Coulomb) given in SAID [21] in Fig. 6. In Table 20 
we show the scattering length a and the effective range r for the nn system in the state 1 5'o. 
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100 200 300 
Ei ab [MeV] 



Figure 6: The calculated nn phase shifts compared to SAID [21] pp phase shift values (without 
Coulomb). 
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